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^ : ABSTRACT 
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^-v , The term "dynamo" means different things to the laboratory fusion plasma 

and astrophysical plasma communities. To alleviate the resulting confusion and 

to facilitate interdisciplinary progress, we pinpoint conceptual differences and 

similarities between laboratory plasma dynamos and astrophysical dynamos. 

^ ! We can divide dynamos into three types: 1. magnetically dominated helical 

dynamos which sustain a large scale magnetic field against resistive decay and 
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O ! drive the magnetic geometry toward the lowest energy state, 2. flow-driven 

Q I helical dynamos which amplify or sustain large scale magnetic fields in an 

r-| I otherwise turbulent flow, and 3. flow-driven nonhelical dynamos which amplify 

5^1 fields on scales at or below the driving turbulence. We discuss how all three 

O I types occur in astrophysics whereas plasma confinement device dynamos are 

'^ I of the first type. Type 3 dynamos requires no magnetic or kinetic helicity of 

^ ' any kind. Focusing on type 1 and 2 dynamos, we show how different limits 

of a unified set of equations for magnetic helicity evolution reveal both types. 

rS I We explicitly describe a steady-state example of a type 1 dynamo, and three 

pi \ examples of type 2 dynamos: (i) closed volume and time dependent; (ii) 

steady-state with open boundaries; (iii) time dependent with open boundaries. 

Key Words: magnetic fields; MHD; Sun: coronal mass ejections; stars: coronae; accretion, 
accretion disks;methods: laboratory; 



1. Introduction 

1.1. General Motivation 

Dynamos describe the amplification and/or sustenance of magnetic fields in electrically 
conducting media such as plasmas and liquid metals. For 50 years the detailed dynamics 
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and mechanisnis of dynamos have been subjects of active research in plasma physics, 
astrophysics, geophysics, and nonhnear dynamics. Historically, much dynamo research 
in these fields has progressed independently. Unsurprisingly, when different communities 
interact, confusion arises when the term "dynamo" is used. This is partly because laboratory 
plasmas are typically magnetically dominated whereas the interiors of astrophysical rotators 
are typically flow dominated. Astrophysicists familiar only with dynamos in flow dominated 
environments wonder what role a dynamo could possibly play in a magnetically dominated 
environment as they are used to thinking of a dynamo as a flow driven amplification of an 
initially weak magnetic field. 

Presently, there is timely motivation to alleviate this confusion and increased 
opportunity for interdisciplinary research. Flow driven dynamos and MHD instabilities are 
being increasingly studied in liquid metals (Gailitis et al.(2002), Gailitis et al.(2003), Ji 
et al.(2001), Noguchi et al.(2002), Peffley et al.(2000), Sisan et al.(2004) to address some 
principles of traditional astrophysical and planetary dynamos, and the magnetically 
dominated dynamos of confinement plasmas are now realized to have direct analogies in 
astrophysical coronae (e.g. Ji et al. 04, Blackman 2005). 

Generally, it is important to clarify what is meant by a dynamo in each context so 
that all research communities can appreciate the common principles and differences. A 
particular unifying question for dynamo theories is to what extent are helical dynamos 
independent of the resistivity (or dissipation). While different answers to this question 
arise in different contexts, we aim to guide the reader toward understanding these different 
answers in a unified framework that draws from recent work on magnetic helicity evolution. 



1.2. Distinguishing Dynamos 

Laboratory plasma dynamos arise in magnetically dominated conditions that have 
been studied in the context of fusion plasma confinement, such as the Reversed Field 
Pinch (RFP) configuration (Bodin & Newton 1990, Ji and Prager 2002). These dynamos 
describe how a large scale magnetic configuration adjusts toward its relaxed state, in the 
presence of external driving away from the relaxed state by a magnetic field aligned electric 
field (or equivalently, external injection of one sign of magnetic helicity) (Ji 1999, Ji and 
Prager 2002, Ortolani & Schnack, Strauss 85, Strauss 1986, Bhattacharjee & Hameiri 
1986, Holmes et al. 1988, Gruzinov & Diamond 1995, Bhattacharjee & Yuan, Bellan 2000). 
The magnetic helicity injection actually drives the system away from the relaxed state, but 
also generates small amplitude fluctuation via kink mode instabilities from large currents, 
or tearing modes from large current gradients. The fluctuations produce a correlation 
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between fluctuating velocity and magnetic fields — the turbulent electromotive force (EMF) 
S = (v X b), which allows the system to evolve back toward a relaxed state. The turbulent 
EMF reduces the field aligned current to restore stability by driving a spatial flow of 
magnetic helicity. This enables the magnetic field structure to evolve toward the largest 
helical scale available (subject to boundary conditions), as this is the lowest energy state 
(Taylor 1986). Continuous injection of magnetic helicity typically leads to a quasi-steady 
dynamical equilibrium with (sawtooth type) oscillations as the system is driven away from, 
and then evolves back toward the relaxed state. If the injection is turned off, the fully 
relaxed state can be reached, but the field eventually resistively decays. Via the dynamo, 
the injection therefore also sustains the large scale helical field against decay. 

In short, laboratory plasma dynamos involve: (1) external injection of one sign of 
magnetic helicity (2) a change the magnetic field structure with conversion of magnetic fiux 
from toroidal to poloidal (or vice versa) geometries, (3) an increase the scale of the field as 
the scale of magnetic helicity increases in the relaxation phase, and (4) sustenance of helical 
field against dissipation. Boundary value studies have focused on all of the above, as well 
as the specific instabilities that drive the small amplitude fiuctuations. Laboratory plasma 
dynamos involve both favorable and unfavorable features for plasma confinement: On the 
one hand they sustain a large scale field in an ordered configuration. But to sustain this 
state, instabilities are required, which produce unwanted dissipation and heat transport. 

The magnetically dominated dynamos of laboratory plasmas can be contrasted to 
the mean field helical dynamo originally applied inside of fiow dominated astrophysical 
rotators (Moffatt 1978, Parker 1979, Krause & Radler 1980, Zeldovich et al. 1983). 
The latter involve an initially weak large scale field which is subsequently amplified via 
strong, large amplitude helical velocity fiuctuations. (In this context, fiuctuations and 
turbulence are used interchangeably.). That laboratory plasma dynamos involve small 
amplitude fiuctuations around a strong large scale field implies that high order correlations 
of fiuctuations can be justifiably ignored. In contrast, the fact that fiuctuations typically 
dominate in a velocity driven dynamo requires a more sophisticated closure for a rigorous 
theory because high order correlations of fiuctuations cannot be straightforwardly neglected. 

An important feature of velocity driven helical dynamos is that the large scale field 
is amplified and sustained on scales significantly larger than the scale of the driving 
turbulence. This contrasts the magnetically driven laboratory plasma dynamo, where the 
fiuctuations vary on time scales short compared to the mean field evolution but can be 
large scale in space. The key similarity between the laboratory plasma dynamos and the 
fiow-driven helical dynamos is that both thrive from a finite S ■ B. This quantity drives 
the large scale helical field relaxation in the laboratory plasma dynamo and amplifies 
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large scale helical fields in the flow-driven case. For the latter, the source of S\\ (where 
II indicates along B) typically, though not exclusively, is the kinetic helicity, (v ■ Vxv), 
a pseudoscalar correlation arising from the interplay between stratified turbulence and 
rotation. For a real system with an outward decreasing density, this can be sustained as 
rising eddies expand and rotate oppositely to the underlying mean rotation to conserve 
angular momentum. Falling eddies rotate in the same direction as the mean rotation. 
Both rising and falling eddies thus statistically contribute the same sign of kinetic helicity 
(or twist) (v ■ Vxv). The northern and southern hemispheres have opposite signs of this 
helicity. Inside astrophysical rotators, the ultimate source of energy for the turbulence is 
typically convection or differential rotation. The magnetic field is not dominant therein but 
the turbulence which facilitates the large scale field growth is necessarily accompanied by 
magnetic fluctuations with energy density comparable to that in the velocity fluctuations. 
Recently, boundary terms, helicity fluxes and anisotropic contributions from mean velocity 
flows have been considered as additional contributions to the electromotive force driving 
the velocity driven helical dynamo, (e.g. Blackman & Field 2000a; Vishniac & Cho 2001; 
Radler et al. 2003, Subramanian & Brandenburg 2004; Brandenburg & Subramanian 2005). 

The growth of small scale magnetic fluctuations in the velocity driven dynamo 
environments of astrophysical rotators also highlights another important concept that 
distinguishes flow driven dynamos from the magnetically driven dynamos in laboratory 
plasmas: laboratory dynamos always involve helicity, whereas both helical and non-helical 
flow-driven dynamos exist. Non-helical dynamos (Kazanstev 1968,Maron et al. 2004, 
Haugen 2003,Haugen & Brandenburg 2004, Schekochihin et al. 2002,Schekochihin et al. 
2004) do not involve a mean turbulent electromotive force, just a turbulent velocity which 
amplifles magnetic energy via random walk line stretching and shear. Non-helical and 
helical flow-driven dynamos also differ in that the former amplifles magnetic energy only 
up to the input driving scale, whereas the latter can amplify flelds on even larger scales, 
as needed to explain observed large scale dynamo cycle periods in astrophysical objects. 
An important complication is that small scale dynamos can operate concurrently with the 
large scale helical dynamo, and the effect of the small scale fleld growth on the large scale 
dynamo has been a subject of considerable research. 

While astrophysical dynamos are typically thought of as flow dominated, it is 
important to note that coronae above astrophysical rotators such as stars and accretion 
disks, are likely magnetically dominated (e.g.Galeev et al. 1979,Schrijver & Zwaan 
2000,Haardt & Maraschi(1993), Field & Rogers 1993). Therefore, astrophysical coronae 
are in fact sites for magnetically dominated dynamos driven by helical fleld injection 
from the astrophysical rotator below (Blackman & Field 2004, Blackman 2005). For all 
extra-terrestrial astrophysical rotators except galaxies, we observe at most the coronal 
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fields, not the interior field. In the sun, stars, and accretion engines there is evidence for 
the presence of large scale coronal magnetic fields. For example, the coronal holes of the 
sun are sites of large scale "open" field lines along which the solar wind propagates (e.g. 
Schrijver & Zwaan 2000). The reversal in sign of these fields indicates that they need to 
be regenerated every 11 years. Analogously, the jets from accretion engines such as young 
stellar objects, active galactic nuclei and even gamma-ray bursts are likely magnetically 
mediated either from fling type models (see e.g. Pudritz 2004) or magnetic towers (e.g. 
Lynden-Bell 2003, Uzdensky & MacFadyen(2006)). The needed large scale coronal fields 
can be produced by the opening up of smaller scale loops from within the rotator below, 
as seen in the sun (Wang & Sheeley 2003). These circumstances reveal that magnetically 
driven helical dynamos, in addition to the flow-driven helical dynamo also play an important 
role in astrophysical magnetic field evolution. Indeed, the most direct analogy of laboratory 
plasma dynamos in astrophysics is the magnetic field relaxation in astrophysical coronae 
subject to helicity injection at the base. The injection drives the system away from the 
relaxed state but instabilities arise, allowing the system to relax. Steady injection wold lead 
to a quasi-steady dynamical equilibrium. 

Recent progress toward understanding helical dynamos has resulted from a combination 
of numerical and analytic work that dynamically incorporates magnetic helicity evolution 
(Pouquet et al. 1976, Kleoorin & Ruzmaikin 1982, Ji 1999, Field & Blackman 2002 
Kleoorin & Rogachevskii 1999, Blackman & Field 2000a, Rogachevskii & Kleeorin 2001, 
Brandenburg 2001, Blackman & Field 2002, Vishniac & Cho 2002, Maron & Blackman 
2002, Blackman & Brandenburg 2002, Kleeorin et al. 2002, Blackman & Brandenburg 2003, 
Blackman 2003, Subramanian & Brandenburg 2004, Brandenburg & Subramanian 2005b). 
However, identification of common principles for laboratory plasma dynamos, helical 
flow-driven dynamos, and magnetically dominated coronal relaxation has been obscured by 
the particular approximations made or by the details of the application. The purpose of 
this paper is to show how a range of different dynamo theories incorporate magnetic helicity 
conservation and how they can be seen to represent different cases of a unified framework. 
We derive the simplest version of the equations representative of a particular type of helical 
dynamo. We purposely do not present complete analyses of the solutions (see Brandenburg 
& Subramanian 2005b for a review); the main goal here is to identify where the equations 
fit into the big picture. 

In section 2 we derive the general magnetic helicity density evolution equations for a 
two-scale system including time dependent and boundary terms. In section 3 we discuss 
how specific cases of these equations correspond to 4 specific dynamo circumstances: (1) a 
steady-state magnetically driven laboratory plasma dynamo, (2) a time-dependent closed 
fiow-driven dynamo, (3) a steady state open fiow driven dynamo and its implications for 
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subsequent magnetically driven coronal magnetic relaxation (4) an open time dependent 
flow-driven dynamo. In section 4 we discuss insights about quenching and steady states 
that emerge from the unified framework. We conclude in section 5. 



2. Unifying Equations for Helical Dynamos 

For completeness we derive the needed magnetic helicity evolution equations and 
consolidate them into a useful form for all subsequent sections of this paper. (See also e.g. 
Bellan 2000). 

We start with the electric field 

E = -V$ - ^tA, (1) 

where $ and A are the scalar and vector potentials. Taking the average (spatial, temporal, 
or ensemble) and denoting averaged quantities by an overbar we obtain 

E = - V$ - dt'K (2) 

Subtracting (2) from (1) gives the equation for the fiuctuating electric field 

e = -V0 - a^a, (3) 

where 4> and a are the fiuctuating scalar and vector potentials. 

Now, using B ■ dtA. = 9i(A ■ B) + E ■ B — V ■ (A x E), where the latter two 
terms result from using Maxwell's equation S^B = — VxE, and the vector identity 
A ■ VxE = E ■ B — V • (A X E), we take the dot product of (1) with B to obtain 

dtiA ■ B) = -2(E ■ B) - V ■ ($B + E x A) = -2(E ■ B) - V • (2$B + A x dtA). (4) 

Eq. (4) describes the time evolution of magnetic helicity density. The same procedure used 
to derive (4) can be used, after dotting (2) and (3) with B and b respectively, to obtain 
equations for the time evolution of the mean magnetic helicity density 



(9i(A ■ B) = -2E ■ B - V ■ (^B + E x A) = -2E ■ B - V ■ (2$B + A x dtA), (5) 
and fiuctuating magnetic helicity density 



dtai-h = -2e ■ b - V ■ (0b + e x a) = -2e • b - V ■ (20b + a x dta). (6) 
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We can eliminate the electric fields from (4-6) by using Ohm's law. Here we consider 
the basic Ohm's law with only a resistive term. For the total electric field, we have 

E=-YxB + r]3. (7) 

Taking the average gives the mean field Ohm's law 

E = -£ - V X B + r^J, (8) 



where £ = v x b is the turbulent electromotive force. Subtracting (8) from (7) gives Ohm's 
law for the fluctuating field 

e = £ — vxb — vxB — Vxb + r^j. (9) 

Plugging (7) into (4), gives 

(9i(A ■ B) = -2r]{J ■ B) - V ■ ($B + E x A) = -2r]{J ■ B) - V ■ (2<I>B + A x dtA). (10) 

Plugging (8) into (5) and (9) into (6) give, respectively, 

dt(A-B) =2£-B-2?7J-B-V-($B + ExA) = 2£ ■B-2r]J -B-V ■ {2WB + Ax dtA) (11) 

and fiuctuating magnetic helicity 

dt'a^= -2e-B-2r]yh-V-((^ + eira) = -28 -B -2r]yh -V ■ {2(^ + a x dtO). (12) 

Dotting (9) with b and averaging reveals the important relation 



£B = eh-r]jh. (13) 

All helical dynamos require a finite £ ■ B. Eqs. (11-13) are the key equations to be used in 
subsequent sections. 



3. Examples of Dynamos from Magnetic Helicity Evolution 

Here we show how different types of dynamos can be understood as limiting cases of 
the equations in the previous section. 



3.1. Steady-State, Magnetically Dominated Laboratory Plasma Dynamo 

In this subsection we work in tlie context of a torus. We take mean quantities to be 
time averages and spatial averages over periodic directions (locally z) and 9, but not over 
radius r (where r = is corresponds to an azimuthal ring at the center of the torus' cross 
section.) While realistic RFPs and Tokamaks are time dependent, to illustrate the relevant 
dynamo simply, we assume mean quantities are steady, and that fluctuations vary on time 
scales much less than the averaging time. For RFPs, sawtooth oscillations and crashes 
occur over millisecond time scales and fluctuations occur in lOOyUs or shorter (e.g. Ji and 
Prager 2002,Ortolani & Schnack). Quantities averaged over > 10ms can be approximated 
as steady. We write the steady-state limit of (12) as 

£ii = ^(v-h-7^rb), (14) 



where h = —0b + ^a x dtS.. Dotting (14) with B, and using (8) gives 



£:B = Vh-?7Jb = ?7JB-EB., (15) 

where, from Eqn. (5), we also have 



E-B = -V-(<I'B + -Ax (9iA). (16) 

The dynamo effect in laboratory plasma configurations such as an RFP emerges when 
a large scale electric field E is externally applied along the initial toroidal magnetic field 
(this represents helicity injection usually via the divergence term, as emphasized below 
(18)). Were there no induced electromotive force £, the measured current term on the 
right hand side of (15) would have to balance the applied large scale electric field along the 
magnetic field E|| in a steady state. But for sufficiently large applied E||, RFP experiments 
reveal (Bodin & Newton 1990, Ji and Prager 2002, Caramana & Baker 1984, Ji et al. 1994) 
that E • B = r] J ■ B > only at a single radius < r = r^ < a, where a is the minor radius 
of the torus and r^. is measured from the toroidal axis. For r < re, E ■ B > ?7J ■ B > and 
for r > Tci r^J ■ B > > E ■ B. Excluding pressure gradient and inertial terms in Ohm's 
law, such measurements imply that S\\ ^ 0. Moreover, since r^J ■ B — E ■ B changes sign 
from negative to positive moving outward through r^ (while J ■ B keeps the same sign), 
S\\ must also change sign from negative to positive across r = r^. Because the third term 
in (14) is often negligible, (14) shows that the divergence of the small scale helicity flux h 
must change sign through r^. In the RFP, the presence of S\\ is sustained by fluctuations 
induced by tearing or kink mode instabilities when the applied |E||| exceeds a critical value. 
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Because averaged quantities of this subsection remain functions of r, the steady 
dynamo just described operates locally in r. Taking the volume integral of (15) and using 
(16) gives 

f£-BdV= fh-dS= f{r]J-B-E-B)dV = f r]3 -BdV + f (^B + ^Ax dt A) ■ dS, (17) 

where we have dropped the third term of (14) as it is typically negligible, and used Gauss' 
theorem to convert the divergence integrals to a surface integrals, keeping in mind that 
for doubly connected topologies we must make sure that h is analytic everywhere. This is 
ensured because our averaged quantities depend only on radius. Accordingly, using vector 
identities, the second surface integral in (17) is 

;$B + i A X 9i A) ■ rfS = f^B-dS--fE-dzfA-rde= /^B-rfS-V,^,, (18) 

^ Jo Zi Jo Jo J ^ 

where Vs is the externally applied voltage drop in the toroidal direction on the surface of 
integration (applied experimentally via gaps) and $5 is the toroidal magnetic flux within 
the surface. On the outer radial surface, there is no normal component of the field and 
the penultimate term would vanish for that surface. The later term of (18) represents the 
helicity injection. This helicity injection is unihelical (i.e. one sign), which is important 
because the magnetically driven dynamo relaxation can be thought of as a process driving 
the injected magnetic helicity the largest scale available subject to the boundary conditions. 

Because of the conducting boundaries, there is no net small scale helicity flow through 
the torus. However, it is instructive to separately consider the inner core (r < Tc) and shell 
(r > Tc) regions (Ji and Prager 2002). The discussion below (16) implies that for r < r^. the 
left side of (17) must be negative. Therefore the second term of (17) must also be negative, 
which from the definition of h, implies an outward flux of positive fluctuating magnetic 
helicity through r^. Analogously, positive fluctuating helicity accumulates into the volume 
defined by r > r^- The helicity flux through Vf. therefore provides the local dynamo. 

Instead of dotting £\\ with B (to get (15)), the dynamo is sometimes expressed by 
dotting (14) with Jm. Using (8) and A = 4^, this simply gives A times (15), so this offers 
nothing new beyond (17) regarding the local nature of the dynamo defined as a sustained 
£\\. However, multiplying (14) by A and taking the volume integral, we obtain 

jSwldV= I KE-BdV = f A{V ■h-7]'yh)dV= f A{7]3 -BdV -E ■B)dV. (19) 



Using 



/ AV ■ hdV = /(Ah) -dS- fh- VAdV, (20) 
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ignoring the third term of (19), and ehminating A in favor of J||, we can re-write (19) as 

I III -IdV = /(Ah)-rfS- fh-VAdV = f {rjjl -E ■3ii)dV, (21) 

where we have used Gauss' theorem to obtain the surface integral, just as described below 
equation (17). Here too the only potentially non-trivial surface terms are radial surface 
terms. If we integrate over the full volume of the plasma (i.e. for all r < a), then the surface 
term in (21) vanishes, as h is measured to vanish at r = a. But (21) then shows that a 
finite VA can produce a global dynamo effect defined by Jtot^W ' J"^^ — ItotV'^iid-V, where 
the subscript "tot" indicates the full volume. The quantity h ■ VA gradient need not be 
finite everywhere for such an effect, only a sufficiently non-zero V ■ h to ensure finite £\\ and 
sufficiently non-zero integral in the third term of (21) are required. For the quasi-local case, 
integration is taken over a sub-range of radii and either of the two middle terms of (21) 
could dominate. Thus / h ■ VAdV ^ emerges as a sufficient but not a necessary condition 
for a quasi-local dynamo. This is consistent with the actual calculation in Bhattacharjee & 
Hameiri 1986), but disagrees with Bellan (2000) as the latter drops the second term of (20). 

For the global dynamo as defined below (21), consider a case in which the magnetic 
helicity associated with the mean field has a locally positive sign in the Coloumb gauge. 

2 

Then A is also positive. If JtotV'^ dV is dominant on the right hand side of (21), then (21) 
implies that the outward fiux of positive small scale helicity (represented by h > 0) is 
anti-parallel to the direction of increasing A overall. Therefore, on average, the dynamo acts 
to homogenize the overall scale of magnetic helicity. This homogenization results because 
the magnetic helicity is injected with one sign, so the dynamo acts to drive the characteristic 
scale of magnetic twist to the largest scale available subject to the boundary conditions. 
Note that for the velocity driven dynamo discussed later in Sec. 3.2, kinetic helicity rather 
than magnetic helicity is injected. There the dynamo acts instead to spectrally segregate 
magnetic helicity of opposite signs while largely preserving a net zero magnetic helicity. In 
general, laboratory plasma dynamos involve the injection and evolution of 1 sign of net 
magnetic helicity (unihelical) whilst kinetic helicity driven dynamos are bihelical. 

The parallel component of the electromotive force can be written S\i = aB. The 
pseudoscalar a for the laboratory case is given from (13) and (15) by 



£B l„:;-ebeibi ,, 

a = -^jT- ~ =2 V ■ h = -^jT- - _o , 22 

B B B B 

where the last similarity follows because the fiuctuations are primarily perpendicular to 
the strong mean fields. The right side of (22) is measured to be larger than any resistive 
contribution (Ji et al. 1994), and the values are consistent with dynamo models based on 
the principles described here. We note that the above descriptions can be extended to the 
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case when global quantities are time dependent (Ji et al. in preparation), as required for 
more precise applications to the RFP. 



3.2. Time-Dependent, Closed, Flow-Driven Helical Dynamo 

In this subsection, instead of considering the averages to be time averages, we consider 
global spatial averages over a closed or periodic volume such that the surface terms vanish 
but the mean quantities remain time dependent. For the laboratory plasma dynamo 
discussed in Sec 3.1, the electromotive force arises via magnetic helicity injection through 
a boundary term. Here kinetic helicity is assumed to be injected into the system and the 
boundary terms are not invoked. These are the circumstances considered in recent analytic 
work (Blackman & Field 2002) and numerical simulations (Brandenburg 2001). Another 
important difference between the laboratory plasma dynamo of the previous subsection 
and the flow-driven helical dynamo here is that the former involves weak fluctuations on a 
strong mean field, whereas the latter involves initially weak fields and strong fluctuations. 

Following Blackman & Field (2002), we distinguish the global volume averages from 
the large scale k = 1 quantities, by using brackets to indicate the former and an overbar 
to indicate the latter. For pseudoscalars, we assume the two averages are equal, (e.g. 
a • b = (a ■ b)). In this case, (11) and (12) become 

dt{a-h) = -2(£.B)-2r]{yh) (23) 

and 



dt{A-B) = 2{£-B)-2r]{J-B). (24) 

In the absence of boundary terms, the magnetic helicity is gauge invariant, but 
choosing the Coulomb gauge we can relate the current helicities to the magnetic helicities 
via (J ■ B) = (/c^A ■ B). To complete the set of equations to be solved, we need an equation 
for S. From its definition, 



dtS = dtvxh + vx dth. (25) 

We now need equations for the fluctuating velocity and magnetic field dfh and (9jv. For 
V ■ V = 0, we have from the induction and momentum density equations respectively. 



(9tb = B- Vv-v- VB + Vx(vx b) - Vxv x b + AV^b, (26) 

and 

dtVg = Pgi{B -Vbi + h- VBi -^^■Vvi + v ■ Vt;^ + b ■ V6i - b ■ Vhi) + uV'^Vq + /g, (27) 
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where f is a divergence-free forcing function uncorrelated with b, v is the viscosity, and 
Pgj = {bqi — V^^VgVj) is the projection operator that arises after taking the divergence 
of the incompressible momentum density equation to ehminate the fluctuating pressure 
(magnetic + thermal). Reynolds rules (Radler 1980) allow the interchange of brackets 
and time or spatial derivatives, so the 5th term of (26) and the 4th and 6th terms in the 
parentheses of (27) do not contribute when put into averages and can be ignored. 

The contribution to dtS from the 3rd term in (25) can be derived by direct use of 
(26) in configuration space. Following Blackman & Field 2002, we assume isotropy of the 
resulting velocity and magnetic field correlations for terms linear in B and also retain the 
triple correlations The contribution to dtS from the 2nd term in (25) also contributes terms 
linear in B, and triple correlations. Here the terms linear in B are best derived in Fourier 
space. For this, Gruzinov & Diamond (1995) invoke the Fourier transform of the terms 
linear in B contributing to (9tV x b, supplemented by a linear expansion of the projection 
operator in k\ << k2, where ki is the characteristic wavenumber of the bracketed or mean 
quantities and k2 is the characteristic wavenumber of the fluctuating quantities b and v. 
Collecting all surviving terms, we then have for (25) 



dt£ = |(b ■ Vxb - V ■ Vxv)B - |v2VxB + z/W x b + Av x V^b + T^ + T^, (28) 



where T^ = v x Vx(v x b) and Tj = (ejg„(Pgj(b ■ Vh — v ■ \/vi)bn are the triple 
correlations. Note that the 3rd, 4th, 6th and 8th terms in (28) come from the v x dfh term 
of (25) and the 2nd 5th and 7th terms come from the dfV x b term of (25). For £11 we have 



dtS^i = (9iV X b + V X dth) ■ B/|B| + V X b ■ (9i(B/|B|). (29) 

Substituting (28) into (29) gives 

dtS\\ = 5bV|B| - /3B ■ VxB/|B| - C^,, (30) 

where a = (l/3)(b ■ Vxb — v ■ Vxv), j3 = (l/3)v2, and ( accounts for microphysical 
dissipation terms, the last term of (29), and T^ + T^ 7^ 0. 

The incorporation of the triple correlations via ( was subsequently named the "minimal 
r" closure (Blackman & Field 2003). In reality, ( can be a function of wavenumber when 
spectral models are considered (Kleeorin et al. 1996). However, the simple minimal 
r closure is an improvement over the first order smoothing approximation (Moffatt 
1978, Parker 1979, Vishniac & Cho (2001)) in which triple correlations are ignored, and 
simpler than the eddy damped quasi- normal Markovian closure (Pouquet et al. 1976). 
In general, triple correlations cannot be ignored because flow-driven dynamos involve 
fluctuations which are generally not small compared to the mean magnetic fleld. Moreover, 
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it is these very triple correlation terms which lead to a turbulent cascade, and the presence 
of a turbulent spectrum. The efficacy of the minimal r closure used for passive scalar 
diffusion by Blackman & Field 2003 was verified numerically (Brandenburg et al. 2004), 
and shows reasonable agreement with numerical tests in the MHD regime (Brandenburg & 
Subramanian(2005)). 

The equations to be solved for the closed, helical flow-driven dynamo case here are 
(23), (24), and (29). The large scale field which grows is fully helical, and so the magnetic 
energy and magnetic helicity equations are essentially the same. The resulting dynamo 
is a modern representation of the a^ dynamo (e.g. Moffatt 1978) that incorporates the 
dynamical backreaction of the magnetic field on the kinetic helicity driving the flow and the 
conservation of magnetic helicity. Solutions of these three equations agree with simulations 
as described in (Blackman & Field 2002). The essence of the dynamo growth is as follows: 
the initial system is assumed to be driven with a finite v-Vxv~ (v-Vxv). This grows 
a finite E\\ which then grows large scale magnetic helicity via (23). But due to magnetic 
helicity conservation, as seen in (24), the small scale magnetic helicity a ■ b must then grow 
in opposite sign to that of the large scale. This also grows the small scale current helicity 
j ■ b of the same sign. This in turn quenches a. The system reaches a steady state: If the 
decay of the large scale helicity were to supersede growth, the small scale helicity would also 
deplete, and growth of the large scale field would again begin if v ■ Vx v is steadily driven. 

In the traditional a^ dynamo of the standard texts (Moffatt 1978), the equation for the 
large scale field is solved, with an imposed form of the electromotive force and the linear 
growth equation results. Blackman & Brandenburg 2003 emphasize that magnetic helicity 
is neither conserved in the equations nor the diagrams of the dynamo in the standard texts. 
In the modern version just discussed, the large scale field evolution equation was replaced 
by the large scale magnetic helicity evolution equation and the additional time dependent 
equations for small scale helicity evolution and turbulent electromotive force evolution are 
coupled into the theory dynamically. 

It should be noted that the two-scale analytic approach has been generalized to a 4 
scale approach (Blackman 2003) to assess whether the small scale magnetic helicity tends 
toward the dissipation scale or toward the forcing scale (with the large scale magnetic 
helicity migrating toward the even larger box scale.) The analysis shows that the small 
scale magnetic helicity ffist appears toward the resistive scale but migrates toward the 
forcing scale before the end of the kinematic regime. Numerical simulations of helical 
dynamos in a periodic box (Brandenburg 2001; Maron & Blackman 2002) also show that 
the magnetic helicity in saturation peaks with opposite signs at the forcing scale and box 
scale respectively. This is a bihelical dynamical equilibrium state in which the small and 
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large scale magnetic helicities both migrate to the largest scales available to them. The 
driving kinetic helicity ensures that these two scales are distinct and prevents the small 
scale magnetic helicity from migrating to the box scale. 

The dynamical backreaction approach of this section can be generalized to an a — Q 
dynamo. Then the large scale magnetic helicity evolution equation must be replaced by the 
vector equation for the large scale field itself, but the small scale and turbulent electromotive 
force equations are also included dynamically (Blackman & Brandenburg 2002). 



3.3. Steady State, Open Flow-Driven Helical Dynamo and Implications for 

Coronal Magnetic Relaxation 

Consider now the limit of (11) and (12) in which the time evolution and resistive terms 
are ignored, but the divergence terms are kept. We then have respectively. 



= 2(£:-B)- V- ($B + Ex A) (31) 

and 

= -2(£-B)- V- (0b + ex a) (32) 

Combining these two equations reveals that the fluxes of large and small scale helicity 
through the system boundary are equal and opposite. This has important implications 
for a helical flow-driven dynamo inside an astrophysical rotator: If helical motions were 
to sustain kinetic helicity inside of the rotator, large and small scale magnetic helicities of 
opposite sign grow as discussed in the previous subsection. The existence of a steady-state 
with open boundaries implies that the boundary fluxes of magnetic helicity contribute to 
the respective loss terms in the large and small scale magnetic helicities. The corona would 
be supplied with bihelical structures (Blackman & Field 2000b, Blackman & Brandenburg 
2003, Brandenburg & Blackman(2003)). This is consistent with time-averaged steady 
coronae of the sun (e.g. Schrijver & Zwaan 2000) and AGN accretion disks (e.g.Haardt 
& Maraschi 1993; Field & Rogers 1993). The bihelical nature of the field supplied by the 
dynamo and the sign dependence of the injected helicity on whether surface shear operates 
on a scale larger or smaller than that of a given loop's footpoint separation are reasons why 
extracting the dominant sign of the solar coronal magnetic helicity in each hemisphere of 
the sun has been somewhat elusive (e.g. Demoulin et al 2002). 

The evolution of magnetic structures injected into a corona is conceptually analogous 
to the evolution of a magnetically dominated laboratory plasmas to injection of magnetic 
helicity such as in a Spheromak (e.g. Bellan 2000, Hsu & bellan 2002). Even though 
the corona in the astrophysical case receives injected helicity of both signs, the guiding 
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principles understood from laboratory plasmas are applicable. The experiment of Hsu & 
Bellan 2002 provides a direct analogy to helical loops of flux rising into an astrophysical 
corona from its rotator below. The loops coalesce at the symmetry axis and form a 
magnetic tower. For large enough helicity injection, the tower can break off a Spheromak 
blob from the kink instability. An astrophysical corona can also be modeled as a statistical 
aggregate of magnetically loops and the corona can be thought of as a single dynamical 
entity (Blackman & Field 2004). The helicity flux to the corona in astrophysical rotators 
acts as a seed for subsequent magnetic relaxation therein. The relaxation opens fleld lines 
that the form coronal holes or jets. 



3.4. Time-dependent, Open, Flow-Driven Helical Dynamos 

In general, both the time dependent and the flux terms in equation (11) and (12) should 
be included dynamically. We briefly describe two calculations of flow-driven dynamos which 
incorporate both, using different sets of approximations. In this section we assume the that 
overbars indicate spatial averages. 

In the context of the Galaxy, Shukurov et al. (2006) have solved the mean fleld 
induction equation for B with £\\ determined from setting dt£\\ = in (2). The £\\ involves 
the difference between the kinetic and current helicities which can be related to small scale 
magnetic helicity in the Coulomb gauge. (Shukurov et al. (2006) formally use a gauge 
invariant helicity density, derived by Subramanian & Brandenburg (2006) to replace the 
use of the magnetic helicity density but the key role of the boundary terms is conceptually 
independent of this.) Effectively, Shukurov et al. (2006) therefore solve the induction 
equation for B, (which depends on £\\ and thus (a • b)) and Eq. (12) for (a ■ b). The 
divergence term in (12) can be replaced with one of the form oc V ■ ((a ■ b)V), where 
V = (0, 0, V z) is the mean velocity advecting the small scale helicity out of the volume. 
This mean velocity also appears in the induction equation for B, highlighting that the loss 
terms in the small scale helicity equation also imply advective loss of mean fleld. This 
approach supports the concept (Blackman & Field 2000a) that a flow of small scale helicity 
toward the boundary may help to alleviate the backreaction of the small scale magnetic 
helicity on the kinetic helicity which drives the dynamo in £\\. However, if V z is too large, 
it may carry away too much of the mean fleld which the dynamo is trying to grow in the 
flrst place. In general, more work is needed to calculate V z from flrst principles, and its 
effect on large and small scale flelds. Coupling the dynamo in the rotator to the magnetic 
helicity evolution in the corona above is also of interest. 

A more restrictive time dependent dynamo that includes boundary terms, maintains the 
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time dependence in (11), but implicitly assumes that equation (12) reaches a steady-state 
has also been studied (Vishniac & Cho (2001)). This approach explicitly incorporates 
the role of shear into the helicity flux. Although the approach involves assumptions that 
have now been avoided in more general calculations of helicity fluxes (e.g. Subramanian & 
Brandenburg 2004) (one being the the flrst order smoothing approximation which can be 
avoided by the "minimal tau" closure discussed in section 3.2), the Vishniac & Cho (2001) 
paper identifies how a time dependent fiow-driven dynamo in a Keplerian shear fiow might 
be sustained by a magnetic helicity flux. We choose to outline this approach in more detail 
here in part because the Vishniac & Cho (2001) is less transparent than the more recent 
paper of (Shukurov et al. 2006) but did flrst explicitly use the boundary flux to solve a 
time dependent dynamo equation. In consolidating the calculation here, we show how it 
dovetails into the unifled framework of section 2. 

If we ignore the resistive terms, the relevant forms of (11) and (12) are 



(9t(A ■ B) = 25 ■ B - V ■ (^B + E X A) (33) 

and 

= -is ■ B - V ■ (0b + elTS) (34) 

Using (34), S\\ can be directly written in terms of the small scale helicity flux as 

B_ .^ B 



S\\ = -—V ■ (0b + e X a) = -—V ■ (-V0 + e) x a. (35) 

Vishniac & Cho (2001) then use (35) in the equation for the mean magnetic fleld applied to 
an accretion disk whose mean quantities are axisymmetric. The mean fleld equation is 

dtB = VxS + Vx(VxB) + XV^B. (36) 

Solving (36) requires use of (35). Vishniac & Cho (2001) invoke a correlation time Tc such 
that a ~ — (e + V(f))rc, (note: they deflne e^f = — e and work with e^j). This reduces the 
last term of (35) to 

2-^V ■ ^^rV0 = --^V ■ Jh, (37) 

which deflnes the helicity flux Jh- (Vishniac & Cho (2001) are missing a factor of 2). 

To proceed, can be inverted in terms of e by Fourier transform, which gives 

dh 



- — eijkej{yi)dkdiei{y. + v), (3^ 



where the spatial derivatives operate on r and an assumption of isotropy in wavenumber 
(not explicitly stated in Vishniac & Cho (2001) was used to obtain this). The notoriously 



-17- 



troubling but useful first order smootliing approximation, wliere terms nonlinear in 
fluctuating quantities are ignored (see section 3.2 above) was then used for e to obtain 
e ~ —V X B — V X b. Ignoring V here in the equation for e (not necessarily justified) and 
assuming r << x Eq. (38) becomes, after some algebra 



JH,i ~ 2/ Tcej{x)dkdiei{y.) = 21 Tceits{BnBtVi{:s.)dndiVs{x.) - BiBtVk{x.)dkdiVs{x.)), (39) 

where / is a spatial correlation scale. Assuming incompressible flow for the fluctuations, and 
assuming that total spatial derivatives of the velocity correlations are small, the third term 
on the right of (39) vanishes and integrating the last term by parts gives 



J^.i ~2/VcB- Vt;iB- Vxv. (40) 

This current is then used in (36) to allow growth of B. Subramanian & Brandenburg (2004) 
and Brandenburg & Subramanian (2005) show that this is one of a number of current 
terms that emerge in a more general calculation which avoids the first order smoothing 
approximation (see section 3.2) with additional fiuxes arising when V is included in e. 
Nevertheless, the importance of Eq. (40) is that it, in principle, allows dynamo growth of 
the large scale field to be driven entirely by the small scale helicity fiux without any kinetic 
helicity. Determining whether this works in practice needs further work. 

The role of the Vishniac-Cho fiux has been investigated numerically in a few 
experiments with somewhat mixed results: Brandenburg & Subramanian (2005a) found 
that the flux can sustain the field growth in the absence of kinetic helicity but only if the 
field already exceeds ~ 70% of the equipartition value of the turbulent field. On the other 
hand, when kinetic helicity is present, the Vishniac-Cho fiux has been numerically shown 
(Brandenbug 2005) to alleviate catastrophic quenching from small scale helicity build up by 
allowing ejection of small scale helicity through the boundary, consistent with conceptual 
suggestions of the role of the boundary terms (e.g. Blackman & Field 2000a). 

The role of the boundary fiux may be particularly important when a rapid, unquenched, 
cycle period is involved such as in the Sun (Blackman & Brandenburg 2003). However, 
as also emphasized in section 3.3, large scale helicity fiux likely accompanies any small 
scale helicity fiux. Significant loss of the large scale field would imply removal of the large 
scale field that the dynamo is invoked to generate inside the rotator thereby lowering 
its maximum value inside the rotator. Care in identifying the relative amount of large 
and small scale helicity fiux is warranted. The same issue also arises in the Galactic field 
calculation of Shukurov et al. (2006). 

A fiux driven helical dynamo is conceptually distinct from the closed volume kinetic 
helicity driven dynamo discussed in section 3.2. The fiux driven dynamo is more analogous 
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to the laboratory plasma dynamo discussed in section 3.1 in the steady state approximation, 
but with a key difference: The source of the helicity flux for the velocity driven case depend 
on fluctuations driven by, and coupled to, the velocity shear (e.g. Balbus & Hawley 1998), 
not fluctuations driven by current driven magnetically dominated instabilities. 



4. Insights from the Unified Framework 

To exemplify the benefit of the unified approach to helical dynamos, we consider 
insights that emerge from synthesizing the cases discussed in the previous sections. In 
particular, we consider what can be learned about saturation, dissipation vs. boundary 
terms, and dynamical equilibria vs relaxed states. As mentioned in the Sec. 1, nonlinear 
saturation and the extent to which helical dynamos depend on dissipation has been a long 
standing topic of study, particularly for flow driven helical dynamos. The magnetically 
driven example of section 3.1 can actually provide insight on the flow driven cases. 

First, the fact that the large scale magnetic fleld dominates the fluctuating fleld in 
section 3.1 highlights that there is no universal requirement that the mean fleld energy in 
the steady state be less than that of the fluctuating magnetic fleld for a working helical 
dynamo. This highlights that the Zeldovich relations (Zeldovich 1956; Zeldovich et al. 
1983) relating the mean and fluctuating flelds when a uniform fleld is imposed and the 
small scale fleld is amplifled non-helically, were never derived as a constraint on helical 
dynamos (Blackman & Field 2005). Second, the laboratory plasma experiments show that 
the resistive terms are small in the RFP dynamos. Eq. (22) reveals that the consistent 
explanation for the sustenance of the mean magnetic fleld via £\\ in an RFP dynamo is 
the flux of small scale magnetic helicity. In a steady-state, the fleld is neither growing nor 
decaying but observed RFP sawtooth oscillations occur on a faster-than-resistive time scale 
and can be accommodated by (22) when augmented by the term ^S^A ■ B from the left 
side of (11). As discussed in the previous section, flux terms can allow the dynamo to incur 
cycles or oscillations on time scales not limited by resistivity. Note the distinction between 
section 3.1 and section 3.2. In the latter, we consider a closed volume time-dependent, 
flow driven helical dynamo. As described therein, the growth starts out fast and formally 
independent of the magnetic dissipation rate. Then as the small scale scale magnetic helicity 
builds up, the dynamo slows to become resistively limited and a steady-state balance occurs 
between the electromotive force and the resistive terms. But, unlike the case of section 3.1, 
the helicity flux terms vanish in section 3.2. Therefore the only terms that can balance the 
electromotive force in the steady state are the resistive terms for section 3.2. 

The role of the helicity flux for a flow driven dynamo is further exemplifled in the cases 
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of section 3.3 and 3.4. When the averages are taken locally, the helicity flux terms for a 
flow driven dynamo can dominate resistive terms. In an open astrophysical object like a 
star or accretion disk the flux terms send magnetic helicity and magnetic energy to the 
corona. Indeed for either the time dependent case of section 3.4 or the steady state of 3.3, 
the electromotive force is balanced by a combination of time evolution of magnetic helicity 
and boundary terms. In both of these cases the dissipation terms are negligible. 

Finally, a further comment on the relation between dynamical equilibria and fully 
relaxed states. When kinetic helicity is steadily injected for a velocity driven dynamo or 
when magnetic helicity is steadily injected for a magnetically driven dynamo, an equilibrium 
or a quasi-steady state can be reached. The latter is exemplified by sawtooth oscillations 
in an RFP dynamo or as the 22 year solar dynamo cycle . The driving or injection take 
the system away from the relaxed state, but the dynamo, fueled by induced instabilities, 
continually competes to take the system back toward the relaxed state. An important 
difference between injecting one sign of magnetic helicity vs. one sign of kinetic helicity is 
that the latter leads to a spectral or spatial segregation of helicities of opposite sign, (the 
bihelical dynamical equilibrium for the o? dynamo discussed in section 3.2). When magnetic 
helicity of one sign is injected, that magnetic helicity seeks the largest scale available to 
it. When the driving is turned off and the resistivity small, the system can fully relax, 
but eventually the field decays. The case of a forced vs. decaying magnetically dominated 
turbulent dynamo was studied in Blackman & Field (2004). Numerical simulations of 
the relaxation of helical MHD turbulence with initially equal kinetic and magnetic energy 
densities were studied in (Christensson et al.(2001)). These works show the tendency of 
a single sign of magnetic helicity to migrate toward the largest scale for magnetic helicity 
injected with a single sign. 



5. Summary and Conclusion 

We have used the equations for magnetic helicity evolution as a unifying framework 
for helical dynamos and we have discussed both magnetically driven laboratory plasma 
dynamos and flow-driven astrophysical dynamos within this framework. We summarize the 
common principles and distinguishing features of these dynamos here. 

Laboratory plasma helical dynamos typically involve a magnetically dominated initial 
state with a dominant mean magnetic toroidal magnetic field. When an external toroidal 
electric field is applied, a current is driven along the field which injects magnetic helicity 
of one sign into the system and generates a poloidal field. For sufficiently strong applied 
electric fields, the system is driven sufficiently far away from its relaxed state that helical 
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tearing or kink mode instabilities occur. The resulting fluctuations produce a turbulent 
electromotive force that drives the system back toward the relaxed state. The relaxed state 
for a such a unihelical dynamo is the state in which the magnetic helicity is at the largest 
scale possible, subject to boundary conditions. When the unihelical helicity injection 
is externally sustained for a real system, a dynamical equilibrium with oscillations can 
incur as the system evolves toward and away from the relaxed state. The time-averaged 
electromotive force is maintained by a spatial (radial) flux of small scale magnetic helicity 
within the plasma. The injection of helicity is balanced by the dynamo relaxation so the 
dynamo sustains the fleld configuration against decay. 

Like the laboratory plasma dynamos, the flow-driven helical dynamos require a 
turbulent electromotive force to grow or sustain magnetic helicity at large scales. These 
dynamos are often invoked as an explanation for the large scale flelds of astrophysical 
rotators. Unlike the laboratory plasma dynamos, for the canonical flow-driven helical 
dynamo, the initial mean fleld is weak and the velocity fluctuations are strong. For the 
simplest time-dependent case in a closed volume, the electromotive force is proportional 
to the difference between the current helicity and the kinetic helicity densities. The latter 
initially dominates and this drives growth of the large scale magnetic helicity by sending 
one sign of the magnetic helicity to large scales and the other sign to scales at or below the 
scale of the dominant velocity fluctuations. Here, unlike the laboratory plasma dynamo, 
no magnetic helicity is injected and so the dynamo acts to segregate magnetic helicity of 
opposite signs spatially or spectrally. The build up of the small scale magnetic helicity 
also grows the small scale current helicity which offsets the kinetic helicity contributions 
to the electromotive force and quenches the dynamo into a steady state. In the absence of 
boundary terms, the steady state is one in which growth is balanced by resistive dissipation. 
When boundary flux terms are allowed, the electromotive force can be sustained by a 
magnetic helicity flux, just as in the laboratory case. Such a helicity flux can arise from 
a Keplerian velocity shear and stratiflcation rather than the current driven instabilities of 
laboratory plasma dynamos. 

In the astrophysical case, when the resistive terms in the magnetic helicity evolution 
equation are small, equal in magnitude (but opposite in sign) small and large scale fluxes 
of magnetic helicity are injected to coronae in the steady state. For a time dependent 
situation, open boundaries can, in principle, allow the dynamo to overcome any long term 
resistive saturation by ejecting the offending small scale helicity from a volume of interest. 
However, such ejection will undoubtedly involve ejection of large scale fleld as well, which 
can reduce the steady state saturation value of the large scale fleld inside the rotator volume 
compared to the asymptotic saturation of the closed case. 
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In the magnetically dominated corona, energy associated with this helicity can be 
extracted into high energy particles. Astrophysical coronae, with their underlying rotators 
acting as magnetic helicity injectors, are the astrophysical circumstance most directly 
analogous to the magnetically driven dynamo physics of laboratory plasma devices such as 
Spheromaks and RFPs. In particular, a single loop with footpoints being sheared or twisted 
provides a direct analogy to the precursor to Spheromak formation. Two subtle differences 
that arise when considering the analogy between coronae and laboratory plasmas more 
carefully are that (1) the corona is composed of many injection sites each analogous to a 
Spheromak, and that (2) the dynamo in the rotator below injects both signs of magnetic 
helicity into the corona. Nevertheless, relaxation of magnetic fields in astrophysical coronae 
can produce the very largest scale fields associated with coronal holes and jets and the 
principles of laboratory plasma dynamos are applicable. 

Our goal has been to provide a base that anchors common principles and differences 
between laboratory and astrophysical helical dynamos to foster further cross-disciplinary 
work. We have avoided a detailed exposition about each type of dynamo in specific 
laboratory and astrophysical systems in order to focus on the basic concepts. 
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